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Structural and electronic properties of CaCu3Ti40i2 have been calculated using density-functional 
theory within the local spin-density approximation. After an analysis of structural stability, zone- 
center optical phonon frequencies are evaluated using the frozen-phonon method, and mode effective 
charges are determined from computed Berry-phase polarizations. Excellent agreement between 
calculated and measured phonon frequencies is obtained; calculated mode effective charges are in 
poorer agreement with experiment, although they are of the correct order of magnitude; and the 
lattice contribution to the static dielectric constant is calculated to be ~40. On the basis of these 
results, various mechanisms are considered for the enormous dielectric response reported in recent 
experiments. No direct evidence is found for intrinsic lattice or electronic mechanisms, suggesting 
that increased attention should be given to extrinsic effects. 



PACS: 77.84.-S, 63.20.-e, 71.20.-b, 77.22.-d 
I. INTRODUCTION 

Recently CaCu3Ti40i2 (CCTO) was discovered to 
possess one of the largest static dielectric constants ever 
measured, reaching nearly eo ^ BO, 000 for single-crystal 
samples at room temperature. Eiu Remarkably, the enor- 
mous static susceptibility is almost constant over a wide 
temperature range ('^0-50010. Moreover, measurements 
of the dynamic susceptibility^ show that e falls from its 
huge static value to a more conventional e ~ 100 with 
increasing frequency, and that this crossover shows a 
Debye-like frequency dependence (e.g., ~ IkHz at 
80 K) with a relaxation rate exhibiting an activated 
dependence on temperature. However, the physical ori- 
gins of the enormous eo and the related Debye relaxation 
process remain mysterious. 

It is tempting to classify a solid exhibiting such a 
prodigious dielectric response as a relaxor or a ferroelec- 
tric, but mounting empirical evidence tends to exclude 
CCTO from either category. The low-temperature static 
dielectric cesponse of perovskite ferroelectrics can sur- 
pass 1000,Q rising rapidly above this value with increas- 
ing temperature and peaking near the ferroelectric phase 
transition temperature. However, this peak in response is 
generally confined to a limited temperature range around 
the transition, whereas that seen in CCTO appears quite 
temperature independent (at fixed low frequency) over 
a much wider range. Furthermore, high-resolution x-ray 
and neutron powder diffraction measurements on CCTO 
reveal a centrosymmetric crystal structure (space group 
/m3) persisting down to 35 K, ruling out a conventional 
ferroelectric phase transition. The Debye- Waller factors 
are also found to be normal over a wide tempjerature 
range, from 35 K to beyond room temperature,lil incon- 
sistent with random local polar displacements of the mag- 



nitude typical in ferroelectrics or their disordered para- 
electric phases. Nanodomains or disorder effects, com- 
mon to relaxor materials, are notably absent: neither 
superstructure peaks nor strong diffuse scattering are 
observed in diffraction experiments. Single-crystal sam- 
ples were describedp|to be mainly twinned, i.e., containing 
nanoscale domainsl differing in the sign of the rotation 
of the TiOg octahedra. The presence of the associated 
domain boundaries may have important consequences for 
the dielectric properties, as we shall discuss further be- 
low. 

The crystal structure of CCTO is body-centered cu- 
bic with four ATiOs perovskite-type formula units per 
primitive cell (where A is either Ca or Cu). The dou- 
bled conventional 40-atom unit cell is shown in Fig. j^. 
The local moments associated with d holes on the Cu" 
cations result in the formation of long-range magnetic 
order; transport, Raman, and neutron diffraction experi- 
ments indicate that CCTO is an antiferromagnetic insu- 
lator having a Neel tenmfti;a|ture Tjv of 25 K (and Weiss 
constant 6^, = -34K).E[0^Eiy Thus, the 40-atom simple- 
cubic cell in Fig. |l| is also the primitive cell of the spin 
structure. Although the extraordinary dielectric behav- 
ior persists far above the temperature scale of the mag- 
netic order, the influence of the strongly-correlated local 
moments, which will likely remain above T/v, is currently 
unknown. 

The purpose of this study is to examine rigorously 
the ground-state properties of this material using first- 
principles calculations in order to identify or exclude pos- 
sible mechanisms underlying the large dielectric response. 
Our calculations within the local spin-density approxi- 
mation can serve as an important first step toward the 
eventual elucidation of this unusual phenomenon. Af- 
ter briefly discussing the technical aspects of our first- 
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FIG. 1. Structure of CaCu3Ti40i2 (CCTO), showing tilted 
oxygen octahedra. White, hght, dark, and black atoms are O, 
Ca, Cu, and Ti respectively. Dashed lines indicate 40-atom 
primitive cell of antiferromagnetic spin structure. 

principles calculations in Sec. II, we provide a detailed 
analysis of the ground-state structural and electronic 
properties in Sec. III. In Sec. IV we present the results 
of frozen-phonon calculations from which we determine 
the frequencies of the zone-center optical modes. As em- 
phasized below, all modes are found to be stable, con- 
sistent with the absence of a structural transition at low 
temperature. Our calculated frequencies agree well with 
those measured by Homes et al.ft although we predict 
one extra mode apparently absent from their data. The 
phonon contribution to the static dielectric constant we 
calculate from first principles does not agree with empir- 
ical observation nearly as well as our computed phonon 
frequencies, but it nonetheless underestimates the mea- 
sured contribution from the JR-active lattice modes seen 
in recent experiments (~80)El by a factor less than three. 
In Sec. V we discuss possible explanations for the ob- 
served enormous static dielectric response and associated 
Debye relaxation behavior, focusing attention on possible 
extrinsic mechanisms as well as considering intrinsic lat- 
tice and electronic mechanisms. We conclude in Sec. VI. 



II. METHODOLOGY 

To investigate the ground-state structural, electronic, 
and lattice-dynamical properties of antifegromagnetic 
CCTO, we employ density functional theoryQ within the 
local spin-density approximation (LSDA)eI via— |the Vi- 
enna ab-initio Simulations Package (VASP);EII13 VASP 
utilizes a plane-pasave basis and Vanderbilt ultrasoft 
pseudopotentialsO Qux pseudopotentials include non- 
linear core correctionslij and, for Ca and Ti, these poten- 
tials treat the highest occupied p shell electrons explicitly 



TABLE I. Comparison of calculated and measured struc- 
tural parameters of CCTO. CCTO has space group Im3 
(point group Th); the Wyckoff positions are Ca(0,0, 0), 
Cu(l/2, 0, 0), Ti(l/4, 1/4, 1/4), 0{x, y, 0). 



Structural parameter 


LSDA 


Exp. (35 K) 


a (a.u.) 


7.290 


7.384 


X 


0.303 


0.303 


y 


0.175 


0.179 



as valence; for Cu we also consider electrons in the 3c? and 
4s shells self-consistently. A 37 Ry plane-wave cut-off re- 
sults in convergence of the total energy to 1 meV/ion. 
The ions are steadily relaxed towards equilibrium until 
the HcUmann-Feynman forces are less than lO^^eV/A. 
When calculating the phonon frequencies, we converge 
the forces to 10~^ eV/A, resulting in precision of less than 
1 cm~^. Brillouin-zone integrations are performed with a 
Gaussian broadening of 0.1 eV during all relaxations. All 
calculations are performed with a 2 x 2 x 2 Monkhorst- 
Pack k-point mesh (equivalent to a 4 x 4 x 4 mesh for 
a single 5-atom perovskite unit cell); doubling the cut- 
off or the size of the k-point mesh does not significantly 
change the total energy or forces. 

III. STRUCTURE 

The structure of CaCu3Ti40i2 was first determined 
frorn-jTieutron powder-diffraction data over twenty years 
ago;li3 these results have also been confirmed by more 
recent neutron-diffraction measurements JH The primitive 
cell is body-centered cubic, with space and point groups 
/m3 and Th, respectively, and contains 20 atoms; this 
space group includes inversion symmetry, precluding the 
possibility of a net spontaneous polarization or ferro- 
electricity in this structure. Because there are three 
Cu^'*' cations in the primitive 20-atom cell, this structure 
cannot easily accommodate an antiferromagnfttic (AFM) 
spin structure. In fact, neutron diffractioifl indicates 
that the AFM state can be described with a doubled 
40-atom simple-cubic unit cell as shown in Fig. 0. In 
this structure, each Cu-Cu nearest-neighbor pair has an- 
tiparallel spins. All calculations are performed using this 
40-atom cell, although it should be kept in mind that 
the atomic coordinates and total charge density have the 
periodicity of the 20-atom bcc cell. 

Beginning with the experimental structural parame- 
ters, the 40-atom cubic cell is first relaxed for a set of 
fixed volumes, and the equilibrium lattice constant is 
then determined by minimizing the total energy as a 
function of volume. Our calculated lattice constant is 
a — 7.290 A, less than the value of 7.384 A measured at 
35 K by slightly more than 1%, as is typical in LDA cal- 
culations. The internal parameters of the structure are 
then relaxed at our calculated lattice constant; their val- 
ues appear in Table |. The stability of the Jm3 structure 
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is also checked by randomly displacing the ions (lowering 
the point symmetry from T/j to Ci) and allowing them to 
relax steadily toward their equilibrium positions. We also 
randomly displaced the ions in the same Cartesian di- 
rection, in this case reducing the point symmetry to C'2v ■ 
Importantly, in both cases the ions returned to their high- 
symmetry (/m3) sites, indicating the absence of unstable 
(i.e., imaginary- frequency) zone-center optical phonons. 
(Here we limit ourselves to lattice distortions having the 
same random displacements in each simple cubic unit 
cell. Thus both zone-center and (111) zone-boundary in- 
stabilities have been investigated for the BCC structure.) 

As is evident from Fig. |^, the /m3 structure of CCTO 
can be obtained from an ideal simple-cubic CaTiOs per- 
ovskite by substituting 3/4 of the Ca ions by Cu in a bcc 
pattern and rotating each TiOg octahedron by a fixed 
angle about one of the four {111} axes. Each Ti ion is 
then fully coordinated with six equidistant oxygen ions, 
calculated to be 1.94 A away. Despite these rotations, the 
Ca ion continues to possess 12 equidistant oxygen neigh- 
bors (calculated distance 2.55 A) as might be expected 
for A cations in the ideal perovskite structure. However, 
the tilted octahedra do result in a radically different local 
environment for each Cu^^ ion: all copper ions are coor- 
dinated by a planar arrangement of four nearest-neighbor 
oxygens (at 1.92 A) and eight oxygens that are consider- 
ably further away. The tilting occurs to relieve tension 
originating from the relatively small ionic radii of both 
A-type cations. Ip_.this context we note that the amund 
states of CaTiOsjIiJ and to a lesser extent SrTi03,E3 also 
exhibit similar octahedral rotations. 



IV. ELECTRONIC AND MAGNETIC 
STRUCTURE 

Our calculation of electronic structure, performed on 
the fully-relaxed geometry, confirms that the ground 
state is an antiferromagnetic (AFM) insulator (at T=0) 
and further illuminates the bonding in CCTO. We find 
the total charge density to be symmetric, and the spin 
density to be antisymmetric, under the fractional lattice 
translation (a/2, a/2, a/2) consistent with the expected 
AFM state. (Here a is the lattice constant.) We plot the 
density of single-particle states (DOS) in Fig. |2[ calcu- 
lated within the LSDA for a single spin channel. Both 
spin channels exhibit identical DOS, as expected for an 
AFM insulator. The bands in the energy range extending 
from —7 to 1 eV are mainly O 2p and Cu 3c? in character. 
The states deeper in these bands, from —7 to— 0.5 eV, are 
found to consist mainly of Cu 3c? orbitals that hybridize 
weakly with O 2p orbitals that point toward their Ti 
neighbors. In the range from — 0.5eV to 1 eV we find 
two narrow bands, each containing just three weakly- 
dispersive states, composed of strongly cr-antibonding 
combinations of Cu 3d orbitals and O 2p orbitals point- 
ing to the Cu from its four near oxygen neighbors (more 
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FIG. 2. Density of states (DOS) for CCTO. Spin-up and 
spin-down DOS are identical, as expected for an antiferromag- 
net. Vertical dashed line indicates valence-band maximum. 

about this to follow below) . The unoccupied conduction- 
band states above 1 eV in Fig. |2| are primarily Ti 3c?-like, 
with mainly t2g states making up the bands in the 1-4 eV 
range and mainly Cg states above 4.5 eV and continuing 
up to about 7.5 eV. 

More details are provided in Fig. H, which affords a 
magnified view of the total DOS near the fundamental 
gap and decomposes the DOS into a site-projected partial 
DOS (PDOS) for the constituent Cu, Ti, and O atoms 
in the crystal. (The PDOS of Ca is quite small in this 
energy range.) By symmetry, the total DOS and the 
Ti PDOS are not spin-polarized, but the spin splitting 
shows up clearly on Cu and O sites where spin-up and 
spin-down densities are indicated by solid and dashed 
curves respectively. Figure || clearly illustrates that the 
fundamental band gap and local magnetic moment arise 
from the splitting of the six weakly dispersive bands of 
Cu(3c?)-0(2p) antibonding orbitals in the — 0.5-1 eV en- 
ergy range. 

Further insight about these states can be gained by 
inspecting Fig. ^, where a contour plot is shown of the 
spin-up and spin-down charge densities associated with 
the lower three occupied valence-band states in this group 
(— 0.5-0 eV); a plot for the three unoccupied conduction- 
band states (0.2-0.7 eV) would look similar but with spins 
reversed. Figure ^ clearly reveals the Cu(3c?)-0(2p) tr- 
antibonding nature of the states, which extend over the 
central cluster composed of a Cu ion and its four close 
oxygen neighbors at 1.92 A (see Sec. HI). This cluster ev- 
idently forms a strongly spin-polarized unit with a ferro- 
magnetic alignment of the magnetic moments on the Cu 
and O atoms comprising it, indicating a significant role 
of the oxygen atoms in the magnetic structure. From 
the spin densities we estimate the magnetic moment of 
each Cu04 complex to be ~0.85/iB, where fiB is the 
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FIG. 3. Total DOS (states/eV-cell-spin) and site-projected 
partial DOS (PDOS) (states/eV-site-spin) for CCTO. Solid 
and dashed curves are spin-up and spin-down densities, re- 
spectively. Vertical dashed line indicates valence-band max- 
imum. The PDOS is calculated using sphere radii of 0.7 A, 
0.7 A, and 1.25 A for Cu, Ti, and O, respectively. The PDOS 
of Ca, not shown, is very small in this energy range. 




FIG. 4. Spin-up (left panel) and spin-down (right panel) 
charge densities associated with the three highest occupied 
valence bands of CCTO plotted with logarithmic contours 
(five intervals per decade) in the z=0 plane. Each density ac- 
tually has mirror symmetry across the central dividing line. 
To facilitate comparison of magnitudes between panels, den- 
sities falling below an arbitrary threshold are shown dashed. 
Atoms are shaded as in Fig. 1. 



Bohr magneton; nearly 85% of this moment is derived 
from electrons in the uppermost three valence bands. 
A closer examination indicates that each of the four O 
ions contribute as much as 0.09 /ig/ion to the total mag- 
netic moment of the CUO4 unit. The development of 
such large magnetic moments on oxygen ions has been 
previously predicted in Li2CU|Q2, a cuprate containing 
one-dimensional Cu02 chains.E3 We also note that our 
calculation compares well with the most recent experi- 
mental value of ^0.8 /is eixtracted from neutron Bragg- 
peak intensities of CCTOjS although less favorably with 
a previous value of ^0.96 //bJ13 Both, however, are sig- 
nificantly larger (by nearly a factor of 2) than moments, 
seen in some undoped cuprates, such as La2Cu04_j,,t3 
which are notably parent compounds of the doped high- 
Tc superconductors. j— . 

We calculate an indirect band gap of 0.19eV££l and 
a minimum direct gap of 0.27 eV. Our computed gap 
thus greatly underestimates the experimental optical 
gap, which must exceed 1.5 eV according to recent data.EHl 
That the LSDA underestimates the gap dramatically in 
this case is not particularly surprising; the localized 3c? 
orbitals centered on the Cu^"*" ions will almost certainly 
result in strong on-site Coulomb interactions, which are 
not adequately treated at the level of the LSDA. In fact, 
we expect these interactions to be sufficiently strong that 
the material is best described as a Mott-Hubbard in- 
sulator. In view of this, some additional details of our 
calculated electronic and magnetic structures should not 
necessarily be taken at face value. For example, we com- 
puted the total energies of the ferromagnetic (FM) and 
unpolarized paramagnetic (PM) states and found them 
to be higher than the energy of the AFM ground state by 
30 and 60meV/Cu ion, respectively. Further, both our 
FM and PM states are metallic. Because of the inade- 
quacies of the LSDA, the temptation to use the AFM- 
FM splitting to estimate the Weiss constant or Neel tem- 
perature should be avoided; both would be unphysically 
large. Moreover, within the conventional picture of a 
Mott-Hubbard insulator we expect the system to remain 
insulating, with little change in the optical gap, as the 
system is switched from an AFM to a FM state, or even 
as the temperature is raised through the Neel tempera- 
ture to create a disordered spin state. 

Clearly, the LSDA has serious limitations for a sys- 
tem of this kind, and therefore it would be of value to 
explore the electronic and magnetic structure of this sys- 
tem using more sophisticated many-body tools. One such 
direction would be along the lines of quasiparticle the- 
ory, perhaps-using the Bcthe-Salpeter extension^ of the 
GW theory.c3 An alternative approach could be along 
the lineS|-af LSDA+Ua or dynamical mean-field theory 
(DMFT)l3. These approaches, which have the capabil- 
ity of taking strong electron-electron interactions into ac- 
count, would presumably give a more realistic description 
of the electronic ground state and its excitations in a ma- 
terial such as CCTO. 
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V. ZONE-CENTER OPTICAL PHONONS 



r = 2Ag © 2Eg ® 4Tg © 2Au ® 2E„ © 12T„ (1) 



In most insulating perovskites such at SrTiOa and 
CaTiOs, the dominant contribution to the static dielec- 
tric constant comes from the lattice susceptibility asso- 
ciated with the IR-active phonon modes. til Thus, we are 
strongly motivated to calculate the frequencies and dy- 
namical mode effective charges of these zone-center op- 
tical phonon modes, in order to test whether the lattice 
susceptibility can explain the enormous observed static 
dielectric constant of CCTO. Moreover, the zone-center 
optical modes are also of interest in view of some seem- 
ingly anomalous behavior observed in the opitijeal IR and 
Raman spectra as a function of temperatureE'cl Only ten 
IR-active modes are observed, despite the fact that the 
symmetry of the crystal requires eleven distinct IR-active 
optical phonons. And, perhaps more crucial, the lowest 
peak, with a frequency of ~122cm~^ at room tempera- 
ture, shifts to slightly lower frequencies, broadens, and 
develops a pronounced shoulder as the temperature de- 
creases. Interestingly, the oscillator strength associated 
with this peak also increases by a factor of 3.6, in appaij^ 
ent violation of the /—sum rule. Raman measurementsH 
reveal not only three sharp phonon lines, but also some 
weaker features and a peculiar continuum extending from 
about 210 to 420 cm^^. It is unclear whether or not these 
IR and Raman anomalies are connected with the huge 
static dielectric response, but it is clearly of interest to 
investigate them as far as possible. 

With these motivations in view, we calculate the fre- 
quencies of all zone-center optical phonon frequencies in 
CCTO, and for IR-active modes we also compute the dy- 
namical effective charges needed to evaluate the lattice 
contribution to the static dielectric constant. 



A. Symmetry analysis 

Since the Hamiltonian remains invariant under a lat- 
tice translation of the 20-atom primitive bcc cell with 
an accompanying spin reversal, the force constant ma- 
trix elements, which are the second derivatives of the 
total energy with respect to atomic displacements only, 
are also unchanged. Thus the force-constant matrix has 
the translational symmetry of the 20-atom cell, and in 
what follows, we use the symmetry- adapted modes de- 
termined from the 20-atom primitive cell (and not the 
40-atom conventional cell) in our analysis of zone-center 
phonons. 

To calculate the phonon contribution to the dielec- 
tric response, we need only determine frequencies and 
mode effective charge of the IR-active phonons at the 
Brillouin-Zone (BZ) center. Altogether there are 60 dis- 
tinct modes at the T point, including the three acoustic 
(zero-frequency) modes. To minimize the numerical bur- 
den, we first decompose these modes into the 6 allowed 
irreducible representations: 



Symmetry-adapted modes transforming as each of 
the irreducible representations given above were con- 
structedoising the SMODES symmetry-analysis software 
package.E3 Among those phonon modes, only the T„ 
modes are infrared (IR) active; the A„ and E„ are silent, 
and the remaining modes are Raman active. 

B. Zone-center phonon frequencies 

For each irreducible representation F, we construct an 
Ny X Ny dynamical matrix from a series of frozen-phonon 
calculations in which the structure is distorted according 
to each of the iVr symmetry- adapted coordinates consis- 
tent with this irreducible representation. For example, 
for the case of the 7Vr=12 IR-active T„ modes (each of 
which is three-fold degenerate), we set up a 12x12 dy- 
namical matrix. In all cases, the forces are calculated on 
the ions in the 40-atom cell and with respect to the cor- 
rect AFM reference ground state within the LSDA. After 
calculating the residual Hellmann-Feynman forces on 
ion i and in Cartesian direction a due to displacements 
Uj of ion j in direction /3, we displace the ions accord- 
ing to the symmetry coordinate for each mode by 0.1% 
of lattice constant and recalculate the forces. The force 
constant matrix 



is then obtained by finite differences by freezing in a small 
amplitude of each symmetry-adapted distortion pattern 
that contributes to the modes of a given symmetry and 
computing the resulting forces that arise. The normal 
modes ux and their frequencies uj\ are then obtained 
through solution of the eigenvalue equation 

*-UA -C^^M-UA. (3) 

Here M^f = (m,/mo) Sij 5a(3 is the dimensionless diago- 
nal mass matrix, where rrii is the mass of atoiB-|Z and toq 
is a reference mass chosen here to be 1 amu,c3 and the 
eigenvectors are normalized according to u^-M-Uj^ = 5^^. 

The results for the IR-active modes are showixin Ta- 
ble The overall agreement with experimental is ex- 
cellent, except for one mode (a; =483 cm"-'^) apparently 
missing in the measurement. By symmetry, there should 
be 11 optical modes (plus 1 acoustic mode) in total, yet 
only 10 optical modes were observed by experiment. It 
is possible that the oscillator strength for this mode is 
immeasurably small; alternatively, its frequency may be 
so close to that of another mode (w =507 cm~^) that the 
pair of modes appears together as a single peak in the op- 
tical data. As we show in the next section, however, we 
compute an appreciable, non-zero oscillator strength for 
this mode, which is mainly composed of displacements of 
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TABLE IL Calculated mode frequencies ujx, effective 
charges ZjJ, and oscillator strengths S\ = ^'^Z'^/uj^ of 
IR-active modes, compared with experimental values ex- 
tracted from Refs. H. 



(cm ^) Zl Sx 



LSDA 


Exper. 


LSDA 


Exper. 


LSDA 


Exper. 


125 


122.3 


0.51 


0.94 


4.1 


14.3 


141 


140.8 


0.65 


1.12 


5.4 


15.9 


160 


160.8 


0.91 


0.84 


8.1 


6.92 


205 


198.9 


0.55 


0.93 


1.8 


5.25 


264 


253.9 


1.79 


1.95 


11.5 


13.8 


314 


307.6 


0.66 


0.52 


1.1 


0.68 


394 


382.1 


0.36 


1.1 


0.2 


1.96 


429 


421.0 


1.83 


1.12 


4.6 


1.72 


483 




1.44 




2.2 




507 


504.2 


0.67 


0.89 


0.4 


0.78 


563 


552.4 


0.73 


0.88 


0.4 


0.62 



TABLE III. Frequencies (in cm~^) of Raman-active and 
silent modes of CCTO, with proposed assignments for the 
three modes observed experimentally in Ref. ^. Irreducible 
representations are given in parentheses. 





LSDA 


Exper. 


Raman modes 


277 (Tg) 






292 (Eg) 






437 (Tg) 






439 (Ag) 


445 




519 (Ag) 


513 




552 (Tg) 






568 (Eg) 


577 




739 (Tg) 




Silent modes 


123 (E„) 
179 (A„) 
472 (A„) 
490 (E„) 





Ti and O ions. Furthermore, the sum of the calculated 
oscillator strengths of the modes at both 483 cm~^ and 
507 cm" ^ would greatly exceed that obtained from the 
observed peak. We return to these inconsistencies below 
in Sec. VI D. 

We also calculate the frequencies of the Raman-active 



modes and silent modes; the results are listed in Table [I] 



In the case of the Raman modes, we also give proposed 
assignments of the experimentally observed modes. In 
the Raman experimento the incidentr-apd scattered light 
were polarized parallel to each other ,EZl in which case the 
symmetry of the Raman-tensor requires the intensity q£ 
all Tg modes to vanish.c3 The experimental spectrumEl 
shows three clear peaks at 445, 513, and 577 cm~^. In 
addition, there is an anomalous temperature-dependent 
continuum extending from about 210 to 420 cm~^, on 
which is superimposed a weak peak at approximately 
290 cm~^. Finally, there is also a weak broad peak at 
approximately 760 cm~^. The three strong peaks are eas- 



ily assigned as in Tabic [II, and the weak low- frequency 



peak is almost certainly associated with the Eg mode 
calculated to lie at 292 cm~^. The feature observed at 
760 cm~^ could be attributed to the Tg mode computed 
at 739 cm" ^ if the assumed polarization alignment was 
imperfect, allowing some leakage of Tg mode intensity 
into the spectrum. We cannot directly explain the contin- 
uum observed in the 210-420cm~^ range, but note that 
it could result from a two-phonon Raman process, possi- 
bly involving pairs of IR-active modes. Overall, we find 
excellent agreement between theory and experiment for 
the frequencies of Raman-active modes. 

To summarize, we have computed the frequency of all 
zone-center modes of CCTO and confirmed that all are 
stable (w^ > 0), consistent with experiments and with 
our structural relaxation. In particular, our computed 
ground state is confirmed to be locally stable against all 
possible ferroelectric lattice distortions. 

To test whether the Im3 structure is truly a global or 



only a local minimum, we also calculated the forces and 
frequencies along the direction defined by some phonon 
modes at different displacement amplitudes, including 
the lowest frequency IR phonon. It is found that these 
phonons remain harmonic for amplitudes up to 4% of 
the lattice constant. We also examined the possibility 
of coupling between the low-frequency 122 cm~^ mode 
(uih) and the Raman modes URaman ((^ =277cm~^ and 
292 cm""'^), which have frequencies roughly twice that of 
the lowest-frequency IR mode, to see if multi-phonon 
processes (i.e., anharmonicity) could explain the anoma- 
lies present in the Raman spectrum and optical response 
measurements. For each of the two lowest-frequency Ra- 
man modes, we calculate the total energy versus dis- 
placement along two separate directions defined by both 
uiR, + URaman and ujR - URaman Combinations of these 
two modes. Our calculations along these four directions 
(two associated with each Raman mode) indicate that 
the total energy is quite quadratic for a wide range of 
amplitudes (up to 4% of the lattice constant), indicat- 
ing that the coupling is too weak among these modes to 
explain the magnitude of the observed anomalies. Our 
calculations do not, however, absolutely rule out the pos- 
sibility of significant coupling between more general com- 
binations of phonons. A more complete stability analysis 
would require the calculation of the full phonon spec- 
trum; we leave this formidable task for a future investi- 
gation. 



C. Effective charges and lattice dielectric response 

We now turn our attention to a calculation of the lat- 
tice contribution to the static dielectric constant, in order 
to test whether the predicted response is of the correct 
order of magnitude to explain the enormous observed sus- 
ceptibility. For this purpose, and also motivated to un- 
derstand better the observed anomalous temperature de- 
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pendence of th e oscillator strength of the 122 cm~^ mode 
(see Sec. VIB ), we calculate the effective charges of the 
IR-active modes and their contribution to the dielectric 
constant. The dielectric function in the frequency range 
near the lattice vibrations can be written as a sum 



(4) 



of electronic and lattice contributions. In most insulating 
cubic perovskite oxides (e.g., ferroelectrics and related 
materials) one has e^o ^ 5 and eo=e(0) ~20-100, so that 
the lattice contribution is expected to dominate. It can 
be written 



eph(w) = ^ 



(5) 



where lux and are respectively the mode frequencies 
and mode dynamical charges, and Oq^ = 4:TTe'^/moV is 
a characteristic frequency having the interpretation of-a 
plasma frequency of a gas of objects of mass mo = l amu,E£l 
charge e, and density (V is the 20-atom primitive 
cell volume). Here, we have used the cubic symmetry of 
CCTO to restrict the sum over A to just one member of 
each of the 11 three-fold degenerate sets of IR-active T„ 
modes, namely the one polarized along z, and the mode 
dynamical charge is then expressed in terms of atomic 
dynamical effective charges Z*,,^ as 



Zt 



(6) 



In practice we compute the mode dynamical charge 
by finite differences as Z'^ = V APz/uq, where AP2 = 
-F2(ucq + uqMx) — ^^(uoq) and uo is a small frozen-ion 
amplitude of the normal mode \i\ added to the equilib- 
rium position of each ion. AP has both ionic and elec- 
tronic contributions. The ionic contribution is obtained 
by summing the product of the displacement of each ion 
with the nominal charge of its rigid core; the electronic 
contribution is computed from firstjjrinciples using the 
Berry-phase theory of polarizationO More specifically, 
the electronic part of the polarization is determined to 
within a constant (an integer multiple of the polariza- 
tion "quantum" eR/T^, where R is a lattice vector) by 
evaluating the phase of the product of overlaps between 
cell-periodic Bloch functions at successive k-points along 
a densely-sampled string of points in k-space. This string 
is chosen so it is parallel to the direction of the desired po- 
larization. For each calculation we use one symmetrized 
string consisting of 4 k-points (corresponding to 8 k- 
points/string in a 5-atom perovskite cell). Increasing the 
number of k-points along the string from 4 to 8 changes 
the effective charges by less than 1%. The electronic con- 
tribution to the polarization is calculated separately for 
each spin channel; the total polarization is then the sum 
of the two spin contributions with the ionic contribution. 

The calculated mode effective charges are listed in 
Table II, where they are compared to values obtained 



from the optical conductivity measurement of Homes et 
alu Surprisingly, though the agreement is excellent be- 
tween our calculated frequencies and those measured at 
room temperature, the theoretical values of the mode ef- 
fective charges are in poorer agreement with the mea- 
sured oscillator strengths. That the experimental val- 
ues refer to room temperature and the computations 
to K may account for some of the discrepancies, al- 
though the large increase measured in the oscillator 
strength of the lowest-frequency mode at low temper- 
ature would compare even less favorably with our cal- 
culations. In this context, we note that the measured 
oscillator strengths are strongly temperature-dependent, 
while the mode frequencies are only weakly so, suggest- 
ing a strong temperature-dependence of the mode effec- 
tive charges. Because the mode effective charges are es- 
sentially electronic quantities, these observations could 
possibly suggest that the electronic structure of this com- 
pound is unusual in some way, and perhaps that this pe- 
culiarity is related to the enormous dielectric response. 
These issues will be discussed further in the next section. 

Having these mode effective charges, we estimate the 
lattice contribution to the static {uj — > 0) dielectric con- 
stant eph(O) = YjX Z'xl^l = Ea 'S'a- Using the eigen- 
modes and their effective charges in Table II, we obtain 
(at zero temperature). Assuming a typical eoo~5- 
10, we expect eo~45-50, to be compared with a sub-IR- 
frequency dielectric constant of ^70 at room temperature 
and ~120 at 10 K. Thus despite some detailed disagree- 
ment between the experiment and theory for the mode 
effective charges, the total contribution of these modes 
to the dielectric constant is of roughly the right order of 
magnitude in this frequency range. 



VI. DISCUSSION AND COMPARISON WITH 
EXPERIMENT 

Although we obtain satisfactory agreement with exper- 
iment for computed structural, electronic, magnetic, and 
lattice-dynamical properties, our calculations do not sug- 
gest any obvious explanation for either the low-frequency 
Debye relaxation or its associated enormous static di- 
electric response. In addition, the observed anoma- 
lies associated with the temperature-dependence of the 
IR and Raman phonons also remain unexplained. The 
strong T-dependence of the oscillator strength of the 
lowest-frequency IR-active mode is especially strange. In 
this section we consider various possible explanations for 
these distinctive features of CCTO and suggest further 
experiments that might shed new light on these puzzling 
phenomena. 
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A. Possible origins of the large dielectric response 
and low-frequency Debye relaxation 

We consider two classes of possible explanations for the 
dielectric response and its frequency dependence, one pri- 
marily intrinsic and one primarily extrinsic. By intrin- 
sic, we mean the response that would be measured in a 
perfectly stoichiometric, defect-free, single-domain crys- 
tal of CCTO. Extrinsic effects, on the other hand, would 
be those associated with defects or other crystal imper- 
fections. 



1. Intrinsic mechanisms 

We first consider the possibility that the origin of the 
dielectric response is actually intrinsic, in which case the 
fact that our calculations do not reproduce the observed 
behavior could be ascribed to some approximation or re- 
striction of our calculation (e.g., the local spin-density 
approximation or the fact that we work at zero temper- 
ature) . 

One class of intrinsic origins for the anomolous re- 
sponse would be related to a lattice instability resulting 
in ferroelectric or relaxor-like behavior. The anomalies 
recently observed in both the IR and Raman spectra in- 
deed raise the possibility of a close connection between 
the lattice modes and the enormous dielectric response. 
Yet, in view of existing experimental evidence and the 
results of our calculations discussed above, any such con- 
nection between the large dielectric response and lattice 
dynamics cannot be described by simple ferroelectricity. 
Additionally, there are several specific objections to a re- 
laxor picture, even from a purely experimental point of 
view. First, relaxors are typically solid solutions (e.g., 
PMN = Pb(Mgi/3Nb2/3)03) in which the random iden- 
tity of atoms on one of the sublattices provides an obvi- 
ous source for disorder. In CCTO the Ca/Cu sublattice 
is well ordered, and hence there is no obvious source for 
relaxor-like disorder. Second, it is very hard to reconcile 
a picture of this kind with the diffraction experiments, 
all of which fail to observe any of the broadening, su- 
perstructure, reduced Bragg intensities, or diffuse back- 
ground that should be associated with such structural 
disorder. Third, either the domains are small enough to 
fluctuate freely in orientation, or they are large enough so 
that reorientation of their polarization can occur only by 
domain-wall motion. In the former case, a temperature 
independent dielectric constant is unlikely. In the latter 
case, the magnitude of the dielectric constant should be 
comparable only to that observed in other ferroelectrics 
containing domains, unless the local polarization is un- 
usually large or the domain-wall mobility monumentally 
facile. 

Moreover, a relaxor picture is difficult to reconcile with 
our theoretical calculations, since it requires the struc- 
ture to be unstable to some kind of local polar lattice 



distortion, and we do not observe such instabilities in our 
calculations. Of course, one could argue that one of the 
IR modes really is soft [oj'^ < 0) even though the LSDA 
calculation indicates otherwise (w^ > 0). However, our 
computed phonon frequencies for the undistorted /to3 
structure are in excellent agreement with the experimen- 
tally observed ones. 

Another possible intrinsic origin of the observed di- 
electric response could be a novel, highly-correlated elec- 
tronic ground state. Given that CCTO is an antifer- 
romagnetic Mott insulator with an electronic structure 
containing relatively localized d states as well as disper- 
sive s and p bands, one possibility would be that its 
anomalous dielectric response could be explaipedJpy the 
formation of electronic ferroelectricity (EFE).c2rE2l As it 
is unlikely that the LSDA calculations would correctly 
describe such a state, our prediction of a conventional 
antiferromagnetic ground state should not be taken as 
definitive evidence against this scenario. Nonetheless it 
can be shown based on symmetry considerations that a 
purely electronic ferroelectric state is not possible. In an 
EFE state, quantities such as the electronic charge den- 
sity and one-particle density matrix would break the Th 
point-group symmetry (and in particular, inversion sym- 
metry). Since the IR- active optical phonons would nec- 
essarily couple to such a symmetry breaking, we expect 
these modes to develop nonzero static mode amplitudes 
in the EFE state, a result that would be evident in the 
diffraction or optical spectra (unless the coupling is un- 
usally small). 

Several more specific objections to other possible elec- 
tronic origins can also be raised. First, an electronic 
instability is most likely in a system with a small gap. 
However, the experiments suggest that the optical gap in 
CCTO exceeds 1.5 eV, as would be consistent with the 
picture of a conventional Mott insulator. It is hard to un- 
derstand how excitons or other low-frequency excitations 
could condense in a system with such a robust gap. Sec- 
ond, regarding our theoretical calculations, it is difficult 
to reconcile our reproduction of the observed phonon fre- 
quencies with the inability of the LSDA to reproduce the 
qualitative nature of the electronic ground state. Third, 
the Debye relaxation time t, which is observed to be ac- 
tivated with an attempt frequency in the MHz to GHz 
range, appears to be too slow to attribute to an electronic 
mechanism. 



2. Extrinsic mechanisms 

In view of the experimental and theoretical difficulties 
with an intrinsic explanation, we observe that many puz- 
zling aspects of the behavior of CCTO become easier to 
understand if one assumes that the large observed static 
dielectric response arises from extrinsic mechanisms, such 
as point, line, or planar defects, or more generally, with 
sample microstructure, morphology, or boundary layers. 
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In this case, the intrinsic static dielectric susceptibihty 
would actually be on the order of ~70-120, as is typ- 
ical Ear conventional dielectrics in the same structural 
class.EJ The large response would then be consistent with 
our theoretical calculations, and with most of the exper- 
imental observations (x-ray and neutron diffraction, IR 
and Raman phonon spectra, etc.) that would effectively 
observe only the conventional dielectric material. In this 
scenario, the observed anomalies in the IR and Raman 
spectra and the huge Debye dielectric response would be 
completely unrelated phenomena. 

Since the dielectric response seems to be too large to 
be explained by point (or line) defects alone, boundary or 
interface effects would be needed. There are two obvious 
candidates for such internal interfaces even in "single- 
crystal" samples: antiphase boundaries separating re- 
gions in which the identity of the Ca sublattice is shifted 
by a primitive-cubic lattice vector, and twin boundaries 
separating domains in which the sense of rotation of the 
oxygen octahedra reverses. Diffraction datacl does not 
rule out the occurrence of either kind of domain wall. 

If such domain boundaries exist, they could possibly 
give rise to a large dielectric response as follows. Sup- 
pose that, away from such boundaries, CCTO actually 
has a weak bulk conductivity resulting from point 
defects of some kind. Moreover, assume that the con- 
ductance of the sample as a whole is eliminated by the 
presence of insulating antiphase or twin boundaries that 
divide the sample into non-percolating, conducting do- 
mains. A static dielectric constant on the order of tins I f 
would then be expected, where eins is the dielectric con- 
stant of the insulating interlayer phase and / is its vol- 
ume fraction. Assuming eins ~ 10^ and / ^ 10^"^, an 
overall static dielectric constant of 10^ is quite plausible. 
Moreover, the observed temperature-independence of this 
macroscopic response is easily explained as long as ejns 
and / do not vary significantly with temperature (in par- 
ticular, a strongly temperature-dependent Ub is allowed). 
Such a model predicts a Debye relaxation behavior, with 
a relaxation time r scaling as , and the observed ac- 
tivated behavior of is then naturally ascribed to an. 
activated behavior of the intradomain conductivity (Tb-Ej 
The observed activation energy of 54 meV could be con- 
sistent with a variety of mechanisms for such an activated 
conductivity. For example, one could assume electronic 
conductivity associated with gap states arising from sto- 
ichiometric variations or dopant impurities, or even ionic 
conductivity associated with mobility of charged defects. 
This model, which is essentially an elaboration of a sug- 
gestion made by Subramanian et al. in a "Note added 
in proof" in th|eir Ref. 1, is worked out in further de- 
tail elsewhereo for various morphologies of the blocking 
boundaries, including the possibility that the blocking 
occurs at the contacts instead of at internal boundaries. 

Perhaps the trickiest aspect of such models is that they 
require the insulating blocking layers to be so uniform 
and robust that they completely prevent conductivity of 
the sample as a whole, despite their small volume frac- 



tion. One would then expect that small variations in 
sample preparation procedures might generate samples 
having enough "holes" in the blocking layers that a mea- 
surable conductivity could be observed. In fact, there are 
hints of such behavior in the experimental work. Ell Fur- 
ther, if the enormous dielectric constant of CCTO does 
in fact originate from the suggested extrinsic mechanism, 
the width of the blocking layers, and thence the size of the 
dielectric response, might be controllable. Accordingly, 
it would be interesting to explore the internal morphoL 
ogy of single-crystal CCTO samples, both theoreticalljS 
and experimentally. In this context, we also note that 
recent studies have uncovered two other new materi- 
als, Ca-doped KT^.^Nb^Og (KTNiCa)!^ and La-doped 
PbTiOs (PLT-A)EI, which also exhibit Debye-like relax- 
ation. However, in KTN:Ca this relaxation occurs at 
temperatures above the ferroelectric transition; in PLT-A 
it occurs at temperatures below the ferroelectric transi- 
tion; and in CCTO it evidently appears in the absence 
of a ferroelectric transition. This would further indicate 
that the relaxation process is not necessarily related to 
ferroelectric or relaxor-like behavior at all. 



3. Suggestions for experimental tests 

Several types of experiments might be useful in test- 
ing whether any of the mechanisms proposed above are 
at work in CCTO. First, it would be interesting to test 
whether, when looking at sample-dependent variations, 
one observes any correlation between the behavior of the 
huge dielectric relaxation (e.g., eo values) and properties 
that are clearly intrinsic (e.g., the oscillator strength of 
the lowest-frequency IR mode) . Such a correlation would 
tend to rule out an extrinsic mechanism and establish a 
connection between the phonon anomalies and the huge 
dielectric response. Second, for samples prepared in such 
a way as to have a measurable DC conductivity, it would 
be very interesting to test whether their conductivity is 
activated in the same way as required to explain the De- 
bye behavior in the macroscopically nonconducting sam- 
ples. Third, C-V measurements might give strong in- 
dications about which type of model is correct. Models 
based on lattice or electronic ferroelectricity would tend 
to predict a saturation of the polarization at a character- 
istic bulk "spontaneous polarization," while the extrinsic 
models predict a more-or-less linear behavior up to some 
kind of breakdown voltage. Nonlinearities would reveal 
the importance of anharmonic effects. 

Some other possible avenues of future experimental in- 
vestigation would be to use electronic probes such as 
angle-resolved photoemission and inverse photoemission 
to study the nature of the valence and conduction band 
edges and to look for unusual quasipartical behavior, or 
to use optical measurements to characterize the dielec- 
tric response in the vicinity of the optical edge. Exper- 
iments that could be sensitive to local structural distor- 



9 



tions (e.g., EXAFS, nuclear electric quadrupole, etc.) or 
to glassy dynamics (NMR, spin echo) would also be worth 
pursuing. An isotope experiment may provide essential 
evidence for or against a lattice mechanism. Experiments 
designed to image the assumed domain-boundary block- 
ing layers (e.g., TEM), and to determine their chemical 
and physical properties, would obviously be of interest. 



B. Experimental IR anomalies 

We return here to a discussion of the anomalous be- 
havior of the low-frequsHcy peak in the IR response ob- 
served by Homes et al.Q This mode, which appears at 
122 cm~^ at room temperature, is observed to broaden 
and shift to lower frequency ('^llScm""'^ at 10 K) as the 
temperature is reduced, where it appears to develop a 
high-frequency shoulder. Most surprisingly, however, its 
oscillator strength is found to increase by almost a fac- 
tor of four as the temperature drops from 300 K to 10 K. 
Moreover, this increase does not seem to be offset by a 
decrease of oscillator strength of other modes; in fact, 
the total oscillator strength summed over all IR modes 
is also observed to rise, growing by ^10% from 300 K to 
10 K. This seems especially puzzling since it appears to 
result in a violation of an oscillator-strength /-sum rule. 

To be precise, it is possible to show that, provided there 
is a good adiabatic separation between the frequency 
scales of lattice and electronic responses, one expects a 
sum rule of the form 



2 f 

TT J 



(7) 



to hold, where the integral runs over the frequency do- 
main of the lattice modes and 



V M,; 



(8) 



has the interpretation of a "plasma frequency" associ- 
ated with the ionic charges and masses. Here i runs over 
atoms in the unit cell of volume V, is the Born ef- 
fective charge (isotropically averagecO) of atom j, and 
the angle brackets indicate a thermal ensemble average. 
This result follows directly from the Kramers-Kronig re- 
lations and the high-frequency lattice dielectric response 
Ree(u;) ~ eoo — i^p/f-^^ expected in the frequency region 
above the lattice modes, where the equations of motion 
are dominated by electric forces and inertia. 

The observed strong temperature-dependence of the 
left-hand side of Eq. (|^) need not violate any physical law. 
Rather, it suggests one of two things. First, it may simply 
be that the atomic Z* values are strongly dependent on 
the local displacements of the atoms, in such a way that 
the thermal expectation value on the right side of Eq. (|^) 
acquires a strong temperature-dependence. Second, if 



some kind of continuum of low-frequency electronic exci- 
tations extends down to the frequency range of the lat- 
tice modes, then the assumptions underlying the deriva- 
tion of Eq. (|^) are not valid. In this case, the electronic 
structure (at fixed ionic coordinates) might be strongly 
temperature-dependent in the temperature range of the 
anomaly. 

Our theoretical calculations do not provide any guid- 
ance in distinguishing between these possibilities. The 
second possibility could be suggestive of some kind of 
electronic mechanism for the enormous dielectric re- 
sponse. However, we note that the temperature and 
frequency ranges of the IR anomalies are rather differ- 
ent from those of the Debye relaxation, so that such a 
connection is not necessarily expected. Perhaps careful 
and precise measurements of e(a;) in the frequency range 
above all phonon modes (i.e., above 750 cm~^) would be 
the best way of testing whether or not the assumed adi- 
abatic separation between lattice and electronic degrees 
of freedom really does occur. Additionally, measuring 
the oscillator strengths as a function of pressure may aid 
in uncovering the nature of the missing mode, as well 
as providing further information about the anomalous 
temperature-dependence of the oscillator strengths of the 
other modes. 



VII. CONCLUSIONS 

We have performed first-principles calculations within 
the local-density approximation on CaCu3Ti40i2 in or- 
der to gain insight into the enormous dielectric response 
reported in recent experiments. Our calculated ground- 
state structural properties and phonon frequencies agree 
quite well with preexisting measurements. The electronic 
structure we calculate appears to underestimate the em- 
pirical optical gap in a manner consistent with previous 
calculations within the LSDA. Our computed effective 
charges are of the correct order of magnitude or better 
and result in a lattice contribution to the dielectric con- 
stant that is roughly consistent with that observed in the 
IR range. 

While this work does not provide a theoretical expla- 
nation for either the Debye relaxation or the unusually 
large and temperature-independent dielectric response 
that is observed, our calculations do appear to limit cer- 
tain intrinsic mechanisms and point toward an extrinsic 
cause. We find CCTO to be stable in a centrosymmetric 
crystal structure with space group /m3, arguing against 
the possibility that CCTO is a conventional ferroelectric 
or relaxor. Our use of the local-density approximation 
constrains our ability to rule out more unconventional, 
purely electronic mechanisms, though we note that there 
is currently no evidence of electronic excitations in the 
optical spectrum. We also add that there is no ex- 
perimental support from existing structural or electri- 
cal measurements for the large dipole moment, whether 
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electronic or displacive in origin, required in each unit 
cell to produce the anomalous response. We are unable 
to provide a quantitative explanation for the changes in 
oscillator strength of the low-frequency IR-active mode 
observed in experiment as a function of temperature, al- 
though, on the basis of existing measurements, we specu- 
late that it may be uncorrelated with the large dielectric 
response. 

An extrinsic mechanism for the static dielectric re- 
sponse and Debye relaxation has many attractive fea- 
tures. In one version of such a mechanism, a weak con- 
ductivity in the interior of domains is barely prevented 
from leading to a macroscopic bulk conductivity by the 
presence of then interface blocking layers. c2l In such a sce- 
nario, the dielectric susceptibility at and below the Debye 
response frequency is controlled by this internal conduc- 
tivity, and is unrelated to lattice distortions, structural 
phase transitions, or novel electronic states of the crys- 
tal. Such a mechanism is then consistent with the ab- 
sence of diffraction evidence for the kinds of symmetry- 
lowering distortions expected if a lattice mechanism were 
at work. It is also consistent with the observation that 
our computed zone-center phonon frequencies are in ex- 
cellent agreement with experiment. By the same token, 
however, such a mechanism cannot explain the observed 
phonon anomalies, which would have to be attributed to 
an unrelated cause. 

In any case, it is of the utmost importance to clar- 
ify which type of mechanism really is responsible for the 
huge dielectric response and the Debye relaxation behav- 
ior in CCTO. For this purpose, the ultimate test must 
be experimental. We have suggested numerous future 
avenues of experimental as well as theoretical investiga- 
tions that could help settle this question. It is our hope 
that the work described here can serve as a firm theoret- 
ical foundation for future investigations of this intriguing 
material. 
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